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Abstract 

Photographic Nuclear Emulsion Detector (PNED) has been in use in nuclear and particle 
physics experiments from the begining, often as the major detector system. However, direct 
measurement of impact parameter in this detector does not seem possible due to some lim- 
itations. This paper discribe a simple yet strong method to estimate the impact parameter 
of events on event-by-event basis. Though this method is develpoed specifically for the pho- 
tographic nuclear emulsion detector, we envision that it should also be applicable to other 
multi-target detector systems. 

Keywords: NUCLEAR REACTION Impact parameter, photoemulsion method, relativis- 
tic nuclear collisions. 

1 Introduction 

A study of relativistic nucleus - nucleus collision is an important tool to produce and investi- 
gate highly dense nuclear matter in the laboratory [1-3]. Based on straightforward geometrical 
considerations, theoretical models predict that the size of the dense nuclear matter zone pro- 
duced in collisions depends strongly on the impact parameter (b) i.e., the transverse distance 
between the center of mass of the projectile and the target nucleus [4]. It is thus very im- 
portant to sort out the collisions according to their centrality. The impact parameter, which 
characterizes the initial state, is not a directly measurable quantity. Thus, it is necessary to 
find out an observable that strongly correlated with it. The simplest observable one can think 
of is the total charged particle multiplicity of an event in case of nuclear emulsion detector. 
We can also try to use the shower (mostlly pions) particle or projectile's proton multiplicity as 
an observable but it does not show a strong correlation with impact parameter due to different 
size of the targets. In emulsion detector; emulsion provides medium to the projectile as well 
as targets. Emulsion mainly composed of H, CNO and Ag(Br). The variation in target size 
could create big confusion in identifcation of events having b = and events having impact 
parameter b > of the bigger size projectile with H and Ag(Br). In this paper, we have devel- 
oped a simple method for estimation of impact parameter on an event-by-event basis, which 
allows one to translate the qualitative estimate of impact parameter from total multiplicity 
into a qualitative one. A quantitative estimate is very convenient in order to present consistent 
results obtained in various experiments. We will also discuss some basic characteristics of the 
interactions with respect to the impact parameter. Some of the charateristics are studied in 
different energy intervals like high, Mid., and low energy. 
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Table 1: Mean Free Path of different projectiles in nuclear emulsion. 



Projectile Energy (A GeV) Mean Free Path (cm.) Ref. 





2.1 


21.80±0.70 


[9] 


l2 C 


2.1 


13.80±0.50 


[9] 


14 N 


2.1 


13.10±0.50 


[9] 


16 


2.0 


12.60±0.50 


[10] 


56 Fe 


1.7 


7.97±0.19 


[11] 


84 Rr 


1.0 


6.76±0.21 


[Present work] 


139 La 


1.2 


5.18±0.30 


[12] 


197 Au 


1.0 


5.60±0.26 


[13] 


238 jj 


1.0 


3.67±0.12 


[14] 


Table 2: The chemical composition of NIKFI BR-2 emulsion [15]. 


Element 




IJJ 12 C 14 N 16Q 


8U Br lU8 Ag 


No of atoms/cc 


x 10 22 


3.150 1.410 0.395 0.956 


1.028 1.028 



2 Experimental Details 

Nuclear emulsion is a detector composed of silver halide crystals immersed in a gelatin matrix 
[5-7] consisting mostly of hydrogen, carbon, nitrogen, oxygen, silver and bromine while a 
small percentage of sulfur and iodine are also present as shown in table 2. In the present 
experiment, we have employed a stack of high sensitive NIKFI BR-2 nuclear emulsion pellicles 
of dimensions 9.8x9.8x0.06 cm 3 , exposed horizontally to u Kr ion at a kinetic energy of around 
1 GeV per nucleon. The exposure has been performed at Gesellschaft fur Schwerionenforschung 
(GSI) Darmstadt, Germany. The events have been examined and analyzed with the help of a 
LEITZ (ERGOLUX) optical microscope having total magnification of 2250X and measuring 
accuracy of 1 /im. In order to obtain an unbiased sample of events, an along-the-track scanning 
technique has been employed [8]. 

The interaction mean free path (A) of 8A Kr in nuclear emulsion has been determined and 
found to be 6.76±0.21 cm. Our collaborators (DGKLMTV Collaboration) [16] found a 
value of mean free path (A) 7.10±0.14 cm. consistent, within the experimental error, with 
our value. We have tabulated the mean free path of different emulsion experiments at similar 
beam energy in table 1. From this table, we may conclude that the mean free path decreases 
with increasing beam mass number at similar energy. The mean free path value obtained in 
our experiment is well fitted in this trend i.e., during event scanning, we picked up all genuine 
events according to our event selection criteria and our criteria of event selection is also right. 
For the present work, we used 1197 events scanned by line scanning method and additional 
162 events were picked up by volume scanning mehtod. The grain density of a singly charged 
particle passing in the same emulsion at extreme relativistic velocity is called the minimum 
grain density (g m i n ). In this experiment, its measured value is equal to 28±1 grains per 100 jum. 
The 8A Kr beam stops within a pellicle. Since, the beam energy decreases as it goes from the 
entrance edge, we have divided each plate in three major energy intervals where the beam has 



energy in the range 0.95 - 0.80 (High Energy), 0.80 - 0.50 (Middle Energy), and below 
0.50 (Low Energy) A GeV, respectively. 

The mean number of fully developed and well separated grains per unit length is called 
the grain density g. It is a measure of the rate of ionization loss. The grain density of 
a track corresponds to a particular specific ionization but its actual value depends on the 
degree of development of the emulsion and the type of the emulsion used. It is therefore, 
necessary to introduce another quantity called normalized grain density which is defined as 
g* = g/g m in- Here g is the observed grain density. All charged Secondaries emitted or produced 
in an interactions are classified in accordance with their ionization, range and velocity into the 
following categories: 

(a) Shower tracks (A,): These are freshly created newly produced charged particles with 
g* < 1.4. These particles have relative velocity (3 > 0.7. For the case of a proton it means 
energy of E p > 400 MeV. They are mostly fast pions with a small admixture of kaons and of 
released protons from the projectile which have undergone an interaction. These conditions 
ensure that showers are filtered from the fragments and knockout protons of the target. 

(b) Grey tracks (N g ): Particles having ionization in the interval 1.4 < g* < 6.0 and 
range > 3 mm are defined as greys. These particles having relative velocity 0.3 < f3 < 0.7. 
They are generally knocked out protons of targets having energy 30 < E p < 400 MeV but also 
admixture of deutrons, tritons and some slow mesons. 

(c) Black tracks (A;,): Particles having range < 3 mm from interaction vertex from which 
they originated and g* > 6.0. This corresponds to a relative velocity (3 < 0.3 and a proton 
with energy E p < 30 MeV. Most of these are produced owing to evaporation of residual target 
nucleus. 

The heavily ionizing charged particles (A^ = N g + jVj, ) are parts of the target nucleus 
and are also called target fragments. 

(d) Projectile Fragments (A/): These are the spectator parts of the projectile nucleus 
with charge Z > 1 having velocity close to the beam velocity. The ionization of projectile 
fragments (PFs) is nearly constant over a few mm and emitted within a highly collimated 
forward narrow cone whose size depends upon the available beam energy. 

The forward angle is the angle whose tangent is the ration between the average transverse 
momentum of the projectile fragments to the longitudinal momentum (pi) of the beam. Taking 
Pl as the beam momentum itself, i.e., 9 F = tan~ 1 (p t /pi J ) =~ 9° in this experiment. The PFs 
are further classified into three categories as follows: 

(i) Heavy Projectile Fragments (A/): PF's with charge Z > 3. 

(ii) Alpha Projectile Fragments (A a ): PF's having charge Z=2. 

(iii) Singly charged relativistic Projectile Fragments (Nf z=1 ). 

Since these PF's have velocities nearly equal to the initial beam velocity, their specific 
ionization may be used directly to estimate their charge. 

The total multiplicity of the secondary charged particle (A^orM) is taken as the 
sum of all charged particles emitted or produced in an interactions (N ch = N s + N h + A/ + 
N a + N fz=1 ). 
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Figure 1: (a) Normalized multiplicity distribution of iVj, at ~ 1 GeV per nucleon. Dotted line 
separate the events from the admixture of CNO target and peripheral collisions with Ag(Br) 
targets, (b) N h distribution of events having N h > 8. Solid line is a double Gaussian fit to 
separate Ag and Br target events. 



3 Method of Target Identification 

The exact target identification in an emulsion experiment is not possible as the medium is 
composed of various elements as mentioned in table 2. However, we can divide the major 
constituent elements into three broad target groups such as H (light), CNO (medium) and 
AgBr (heavy) with high accuracy. There are a lot of other ways of statistical separation [17- 
20] , which roughly give the probability of interactions with different targets. It is well known 
that the number of heavy particles, N h is a good tool for target identification. Since we are 
interested in the separation of targets on event by event basis, we have employed short-range 
track distribution to identify targets for events with low Nh value. In view of the distribution of 
heavily ionizing charge particles for all set of events as shown in fig. 1(a), we have attempted 
the separation of targets using the following criteria: 

H target events : Nh — 0; Nh — 1 but not falling in any of the below categories. 

CNO target events : 2 < N h < 8 and no track with range > 10 /jm. 

Ag(Br) target events : N h > 8; N h < 8 and at least one track with range > 10 /iin and 
no track with 10 < range < 50 /im. 

As a result, we have obtained the percentage for the occurence of the three different target 
group events as H : 12%, CNO: 48% and Ag(Br): 40%. The relevant data is summarised in 
table 3. Table 3 gives the results using the above criteria along with the results of other 
similar efforts [21-28]. It may be seen from the table that the probability of events due to 
Ag(Br) nuclei increases slowly with increasing projectile mass at similar energy. It also shows 
that the method of target separation is almost correct. 

We are, first-ever, separating the Ag and Br nuclei in emulsion detector. For this, we 
examine only Ag(Br) type interaction. Their frequency distribution is shown in fig. 1(b) and 
has been fitted by the double Gaussian functions. We have found that the contributions of 



Table 3: Percentage of interactions with different target groups. 



Interactions 


Energy (A GeV) 


H 


CNO 


Ag(Br) 


Ref. 


l P 


2.5 


18.00 


49.50 


32.50 


[21] 


A He 


3.7 


21.03 


40.42 


38.55 


[22] 




3.7 


21.29 


30.87 


47.84 


[23] 


22 Ne 


3.7 


12.94 


32.59 


54.47 


[24] 


28 Si 


3.7 


15.29 


33.79 


50.92 


[25] 


m Ar 


1.8 


17.80 


34.60 


47.50 


[26] 


56 Fe 


3.7 


23.13 


22.64 


54.23 


[27] 


84 Rr 


1.0 


12.10 


47.60 


40.40 


[Present work] 


197 Au 


8.7. 


19.00 


36.00 


45.00 


[28] 




Figure 2: Reference distribution of b and total charge particle multiplicity for different target 
groups and for different energy intervals. Solid line is just to guide the eye. 

Ag nuclei and of Br nuclei are 44% and 56% , respectively. The same numbers for Ag and Br 
nuclei are reported by B. Jakobsson et al [29], who employed the distribution of the sum of 
charges in events where all targets have been measured directly from an exposure at GANIL. 
From now, we will proceed for development of new method for impact parameter estimation 
on event by event basis only for separate target group to avoid the mathematical complications 
and making the method simpler. 

4 Method of Impact parameter (b) Estimation 

Shower particle multiplicity have strong correlation with impact parameter in the fix target 
mass experiments. Due to the invariant target mass, shower particle multiplicity will not show 
strong correlation with impact parameter. Because H target, with fix mass beam, having 
nearly zero impact parameter may show the similar shower multiplicity as CNO or Ag(Br) 
target, with fix beam, having larger impact parameters. That's the reason we are interested in 
using total charged particle multiplicity (N c h) to estimate the impact parameter and will show 
strong correlation between them. In some collider experiments [30] they found that the total 



particle multiplicity (N ch ) is strongly correlated with the impact parameter b. More precisely, 
its mean value decreases monotonically as a function of b. 

In this method the basic assumptions are as following: We treated interactions of each 
target group with projectile separately to make the method simpler. The total cross section 
is purely geometrical and there is a strong correlation between b and N c h- The targets are 
randomly distributed in the nuclear emulsion detector and the probability of interactions is 
totally random and the collisions geometry is also random. 

Therefore, we used random number generator [31] to generate random numbers in be- 
tween and 1 for N c h and b. For making the random number values upto the order of total 
multiplicity {N c h), we choose the maximum and minimum limit according to our real ex- 
perimental values of maximum and minimum multiplicities of total charged particles in each 
target group and each target group multiplied by a suitable factor of numbers. For making 
the random number values upto the order of impact parameter according to the Ap and At, 
we used real radius values of projectile and target nuclei and took the maximum and mini- 
mum as a combination and difference of radii of both nuclei, for b = Rp + Rt = b max and 
b = \R P ~ R T \ = 6 m j„(~ 0), respectively and made the order of magnitude according to the 
order of real experimental values. We assume that minimum multiplicity value belong to the 
b max and maximum multiplicity belongs to the b m i n and therefore we set two extreme points 
on the b versus multiplicity distribution and distributed rest of the values according to the 
monotonic relation between b and N c h such as b = N ch ^ max ) - N c h for different target groups 
and for different energy intervals. In fig. 2, we are showing relation between generated impact 
parameter and generated multiplicity of charged particles. We know the real event total charge 
particle multiplicity N c h and match this real multiplicity with b with the help of fig. 2 and 
can easily estimate the impact parameter corresponding to the real event. 

5 Experimental Results 

To check the authenticity of the developed method, we have checked few charecteristic pa- 
rameters with respect to impact parameter. Most of them are also checked in different energy 
intervals as well as different target groups. In this sequence, first we checked total charge 
particls multiplicity with respect to impact parameter for different target groups in different 
energy intervals as shown in fig. 3 (a), (b) and (c), respectively. The error bars shown 
on the data points are purely statistical and different type of lines are the best fit of data in 
different energy intervals. From fig. 3, we may infer that in each target group, the relation 
between multiplicity and impact parameter is monotonic and the slope is variable with change 
in beam energy as well as change in the target mass number. The nature of the distribution is 
similar to the participant - spectator Model prediction i.e., at maximum overlap (b min ) region 
of target and projectile, multiplicity must be maximum and vice - versa. 

We checked the emission of alpha particles in interactions with H - target with respect to 
the impact parameter as shown in fig. 4(a). The distribution shows a polynomia nature. 
The maximum number of emitted alpha particle in an interaction with H - target is 5 to 6 
at maximum overlap. Here the size of H - target is much smaller than the size of the 8A Kr- 
projectile. That's why it is very difficult to select head-on collisions and this can easily mix 
into the central as well as quasi-central events. Emission of the average number of shower 
particle versus b is plotted in fig. 4(b). The condition of this distribution is similar to the 




Figure 3: The average total charge multiplicity is plotted as a function of estimated impact 
parameter in different energy intervals for different target groups (a) H, (b) CNO, and (c) 
Ag(Br). Solid, short dash and long dash lines are the best fit of high energy, mid. energy and 
low energy data set, respectively. 

fig. 4(a). It shows a strong monotonic correlation between < N s > and b. The solid line is 
the best fit of data points and the error bars represent the statistical error. 

The Q (= £/z=i + 2 x S/z=2 + x fz>?) is the total projectile fragments charge flow 
in the forward emission cone. The < Q > value versus b for all energy events in different 
target groups is shown in fig. 5. For making a comparison and guiding the eye, we have 
fitted the distribution with a linear function. The light and heavy target groups show nearly 
similar slopes while medium target group has slightly different slope. We can also see the 

< Q > value change with target mass, specially with light and heavy targets. At b m i n and 
bmax, the difference in < Q > is around double. At b m i n , less than half of the projectile mass 
are converted into some other things (like energy and neutral particle) but at b max nothing is 
going to disappar while few neutrons are converted into protons and make the total charge Q 
value larger than 36. 

In fig. 6, we have plotted < b > versus Q for all high energy events for one combined 
target (H + CNO + Ag and Br). It can be seen that Q is minimum at minimum value of 

< b > and approaches the beam charge (Z = 36) at maximum value of < b >. The presence 
of some events having Q values more than the value of beam charge (36) when < b > has 
a value less than the maximum < b > may be attributed to the conversion of neutrons into 
protons while nuclei collide with each other. All above described behaviour of paramerets are 
consistent to the several theoretical models [32] specially participant spectator model [33]. 

6 Conclusions 

A new method to estimate an Impact parameter on the basis of event-by-event in Photographic 
Nuclear Emulsion detector, was developed and studied in detail of some characteristic param- 
eter. Although this work originates from the goals of handling data from photographic nuclear 
emulsion detector, and can be applicable to the similar kind of detectors or other multi-target 
detector systems where all type of targets are mixed together and only very limited informa- 
tion is in hand. These studies will boost the data analysis for emulsion detector in the light of 
impact parameter. This is a new method that's why we have no any other experimental data 
for comparison but the experimental values obtained in this experiments strongly support the 
theorectical concepts. 




Figure 4: (a) Average number of Helium nuclei emitted in interaction with H - target are 
plotted with respect to impact parameter for all energy. Solid line is the polynomial function 
fit is just to guide the eye. (b) Average number of shower particles emitted in interaction with 
H - target are plotted as a function of impact parameter for all energy range (0.06 - 0.95 A 
GeV). Solid line is the best line fit of data points. 




Figure 5: The average Q values versus impact parameter for high energy interval with different 
target groups. Different lines are best fit of the data points and are guiding to the eye. 




Figure 6: The average impact parameter versus Q for high energy interval. Solid line is the 
best fit of data points and is just to guide the eye. 
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